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A hypothesis on the physical mechanism generating the heartbeat instability in complex (dusty) 
plasmas is presented. It is suggested that the instability occurs due to the periodically repeated 
critical transformation on the boundary between the microparticle-free area (void) and the complex 
plasma. The critical transformation is supposed to be analogous to the formation of the sheath in 
the vicinity of an electrode. The origin of the transformation is the loss of the electrons and ions 
on microparticles surrounding the void. We have shown that this hypothesis is consistent with the 
experimentally measured stability parameter range, with the evolution of the plasma glow intensity 
and microparticle dynamics during the instability, as well as with the observed excitation of the 
heartbeat instability by an intensity-modulated laser beam, tuned to the atomic transition of a 
working gas. 
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I. INTRODUCTION 

Instabilities in electronegative discharges are known 
since a long time [T]. The so-called attachment instabil- 
ities in glow dc discharges [2], as well as in capacitively- 
coupled [3] and inductively-coupled [|1 rf discharges were 
observed. In the latter case the instability was found to 
be accompanied by the E-H transition |5I6) . All these 
instabilities reveal themselves as relatively low-frequency 
(1-100 kHz) self-excited oscillations of electrical and op- 
tical parameters of the discharge. Variation of opti- 
cal emission has a global character, i.e., the emission 
changes in phase in the whole discharge. Therefore a 
global description implying homogeneous distribution of 
all plasma species is sufficient for understanding of such 
instabilities. 

Complex (dusty) plasmas can also be considered as 
plasmas with very heavy negative ions. However, due to 
the macroscopic nature of microparticles, complex plas- 
mas exhibit two features which make them remarkably 
distinct from regular electronegative plasmas: (i) The 
charge of a microparticle is determined by the dynamic 
balance of the ion and electron fluxes on its surface [7], 
which represents therefore a volumetric sink for a plasma, 
(ii) The ion drag force, acting on microparticles, becomes 
important, leading to the formation of so-called voids, 
i.e. microparticle-free regions [HEl]. Such plasmas can 
no longer be treated as homogeneous. 

In the first microgravity experiments with complex 
plasmas, a spontaneous periodic contraction of the void 
boundary was reported |15j . Because of its characteris- 
tic appearance as well as due to very low repetition fre- 
quency (from single contractions to hundreds Hz), this 
phenomenon was termed "heartbeat" instability [TBlll8j . 
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Several years after the discovery, it was shown in ground- 
based experiments that the heartbeat instability can- 
not be reduced to the dynamics of microparticles only 
[inillT]. The global characteristics of the discharge (e.g., 
rf current) were found to be modulated with the low fre- 
quency of the instability. High-speed real time imaging of 
the discharge showed that the contraction of the void is 
preceded by a steep increase of the glow inside it. Physi- 
cal mechanisms leading to such a behavior of the plasma 
are not understood until now. Recent review of exten- 
sive data from the microgravity experiments, conducted 
on board of the International Space Station [T^ , also did 
not bring any new insights into the understanding of this 
phenomenon. 

In this paper we suggest a hypothesis on the physi- 
cal mechanism of the heartbeat instability in complex 
(dusty) plasmas. This hypothesis employs both distinct 
features of complex plasmas - heterogeneity and the pres- 
ence of a volumetric sink of plasma - and does not contra- 
dict to any available experiments. Under this approach 
the driving mechanism of the heartbeat instability is the 
formation of a sheath at the void boundary, leading to 
the abrupt change in the profile of the electrostatic po- 
tential. 



II. EXPERIMENTAL SETUP 

Our experiments were performed in the PK-3-f cham- 
ber |18j . which is a symmetrically driven parallel plate 
capacitively coupled rf reactor (see Fig.[T]). The discharge 
was created by applying 13.56 MHz voltage to the disc- 
shaped electrodes of 6 cm in diameter, separated by a 
3 cm gap. As a working gas we used argon in a pressure 
range oi p — 10 — 50 Pa. 

Melamine formaldehyde plastic spheres with the diam- 
eter 2a = 1.95 /im were injected into the discharge and 
levitated in the vicinity of the bottom electrode. By heat- 
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FIG. 1: Experimental setup for the investigations of the heart- 
beat instabihty. Microparticles are levitated in the discharge 
volume with the help of the thermophoretic force. A 532 nm 
laser is used for the illumination of the microparticles. A 
772 nm tunable laser, modulated by a chopper, provides ad- 
ditional means of exciting the heartbeat instability. Attenu- 
ators of different transparency are used to vary the tunable 
laser beam power. Two photomultiplyers are recording time- 
resolved emission intensity of the plasma. Longpass filter in 
front of them rejects the scattered light of 532 nm illumination 
laser. A CMOS camera records the dynamics of microparti- 
cles. 



ing up the bottom electrode and thus compensating the 
gravitational force by thermophoresis [2U] a significant 
part of the discharge volume was filled with microparti- 
cles, forming a void almost in the center of the chamber. 
The temperature difference between the electrodes was 
19 K. The microparticles were illuminated by a vertical 
sheet of green (532 nm) laser light and observed from 
the side with a CMOS video camera at a frame rate of 
1000 fps. 

At two radial positions - one inside the void and the 
other outside it, on the periphery of the microparticle 
cloud (see Fig. [T]) - the light from the plasma was 
collected by two lenses and by means of optic fibers 
transmitted to the two respective photomultiplier (PMT) 
modules with an amplifier bandwidth of 200 kHz. 

Along with the observation of the sporadically trig- 
gered heartbeat instability, we used a tunable diode laser 
to excite it " on purpose" . The wavelength of the laser was 
set to the center of the Doppler profile of A = 772.38 nm 
spectral line of argon. This light is resonantly absorbed 
by argon atoms in the IS5 metastable state, causing the 
excitation into the 2p7 radiative state. The 2.4 mm di- 
ameter cylindrical beam from the tunable laser with the 
maximal power of 8 mW was modulated by a mechanical 
chopper and then horizontally driven through the cham- 
ber. It was possible to move the beam vertically. The 
laser power entering the chamber was varied by attenuat- 
ing the beam with the neutral density filters. The same 
laser can be used for measuring the number density of 
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FIG. 2: Example of the stability diagram, showing the param- 
eter range where the heartbeat instability is triggered (inside 
the dashed area) for a given number of levitating microparti- 
cles. Circles represent the measured stability boundary which 
can vary significantly between different experimental series. 

Is5 metastable argon atoms by optical absorption in a 
manner similar to that of Ref. [21] 

III. RESULTS 
A. Reproducibility 

A special note should be made regarding the repro- 
ducibility of the results. We were unable to get rid of the 
uncontrolled drift of the instability parameters. Signif- 
icant change in, e.g., frequency of the heartbeat, some- 
times observed within 10 — 15 min interval, while both 
discharge power and pressure were kept constant. There- 
fore, in each series of measurements we tried to keep the 
measurement time as short as possible, and always dou- 
ble checked that the observed effect is mainly due to the 
variation of the control parameter and not due to the 
degradation of the discharge conditions. 

B. Self-excited heartbeat instability 

The presence of the instability was found to be deter- 
mined by a combination of three parameters: pressure, 
discharge power and number of microparticles. An exam- 
ple of a 2D stability map (at fixed number of microparti- 
cles) is shown in Fig. [2] Obviously, there exists a thresh- 
old pressure (47 Pa in the shown example) , above which 
no instability can be excited at any discharge power. At 
lower pressures the instability can be triggered in a cer- 
tain range of discharge powers, which broadens as the 
pressure is reduced. Below 10 Pa the discharge with mi- 
croparticles cannot be sustained. For given pressure and 
power the instability sets in easier when the number of 
microparticles is larger. 
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FIG. 3: Evolution of the integral light intensity measured 
during the self-excited heartbeat instability (a) outside the 
void and (b) inside the void. The curves are obtained by 
averaging over 200 consecutive cycles of the instability. 

In Fig. |3] the evolution of the intensity of the plasma 
emission inside and outside the void is illustrated. To re- 
duce the noise curves are averaged over 200 consecutive 
heartbeat events. We also made sure that two instability 
cycles never occurred within the time span of Fig|3] Fig- 
ure [3] confirms the earlier observations [T71 [TH] . During 
the void-collapse phase of the heartbeat instability the 
glow intensity first rapidly (about 0.2 ms in our case) in- 
creases inside the void and simultaneously decreases out- 
side. Then both intensities exhibit a plateau for several 
ms. This is followed by an overshoot phase, when the 
intensities go above and below the average level outside 
and inside the void, respectively. Relaxation of this over- 
shoot lasts until the next cycle of the instability starts. 
We would like to stress here that the intensities inside 
and outside the void oscillate remarkably out of phase. 



C. Laser excitation of the heartbeat instabihty 

Before performing the laser excitation of the heartbeat 
instability the system was brought to a stable condition. 
This was obtained by decreasing the discharge power 0.1- 
0.2 W below the lower onset threshold (which was ob- 
served between 0.8 and 1.0 W for the laser excitation 
experiments). Then, by shining an intensity- modulated 
light from a laser tuned to 2p7 — >■ IS5 transition of argon 
we were able to observe the heartbeat instability [22] . We 
note, that the laser excitation of the heartbeat instabil- 
ity was performed in a dedicated series of experiments, 
conducted at p = 20 Pa. This explains a significant devi- 
ation of the lower power threshold of the instability from 
that in Fig. [2| 

The excitation of the heartbeat instability by a mod- 
ulated beam of a tunable laser exhibits resonance prop- 
erties. It was always observed in quite a narrow band of 
the beam modulation frequency fias around 40 — 60 Hz. 
Example of a resonance curve is shown in Fig. [4^. We 
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FIG. 4: Dependence of the frequency of the laser-excited 
heartbeat instability, fub, on (a) the frequency of the laser 
modulation fias and (b) the tunable laser power Pias- The 
results are for (a) Pias ~ 3.5 mW and {h)fias = 55 Hz. Laser 
beam crosses the center of the void. 



see that the frequency of the heartbeat, fhb, can be used 
as a measure of "strength" of the instability. The laser- 
excited heartbeat instability becomes stronger with the 
increase of laser power, Pias (Fig. |4|d) . 

The heartbeat instability strength depends also on the 
position of the laser beam f22]. Figure [s] shows the depen- 
dence of the frequency of the heartbeat instability on the 
fraction of the void volume occupied by the laser beam. 
The heartbeat frequency is maximal if the laser hits the 
center of the void, and decreases as the laser moves out 
towards the periphery of the cloud. 

The heartbeat instability was also found to disappear 
if the laser is detuned from the argon atomic transition. 



IV. DISCUSSION 

A. Critical transition at the void boundary 

The importance of the boundary between the void and 
surrounding complex plasma is substantiated by three 
experimental facts: (i) the instability never occurs with- 
out a void, (ii) variations of the intensity of the plasma 
glow are opposite in phase inside and outside the void, 
(iii) laser excitation of the instability is only possible if 
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FIG. 5; (a) Positions of the laser beam with respect to the 
void, at which (b) the dependence of the frequency of the 
laser-excited heartbeat instability, fht, on the fraction of the 
void volume occupied by the laser beam, V/Vv, was investi- 
gated. The void of volume Vv is approximated as an oblate 
spheroid with the large and small semiaxes 6.5 and 2.6 mm 
respectively. Laser beam is a cylinder with the diameter of 
2.4 mm, V is defined as a volume of intersection of the laser 
beam and the void. For the shown example, Pias = 8 mW 
and fias = 50 Hz. Movies of the microparticle cloud with the 
laser beam in the center of the void (V/Vv = 0.12) and for 
the laser beam outside the void {V/Vv = 0) can be viewed in 
Ref. [22]. 

the plasma parameters are modulated inside the void. 

In low-pressure discharges the electron impact ioniza- 
tion is usually counterbalanced by the loss of electrons 
and ions on the electrodes and walls of the vacuum cham- 
ber. As pointed out above, microparticles immersed into 
a gas discharge provide an additional - volumetric - loss 
of charged particles. Loss rate on microparticles is pro- 
portional to their number density Ud- 

Let us consider a boundary between a void and com- 
plex plasma with so small nd, that plasma losses on the 
microparticles are negligible. In this case the electric 
field is simply determined by the ambipolar diffusion of 
charged particles to the electrodes and walls of the vac- 
uum chamber. In an opposite limiting case can be 
considered so high that the presence of microparticles is 
equivalent to the placement of a solid floating electrode 
enclosing the void into the discharge. In this regime a 
sheath should be formed at the void boundary. 

Transition from a smooth potential profile to forma- 
tion of a sheath around the void should therefore occur 
at a finite value of Ud when the ion loss rate on the mi- 
croparticles is so large that the ions at the void boundary 
would have to be accelerated to a critical velocity (which. 
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FIG. 6; Schematic representation of the sheath formation 
near the void boundary. Abrupt steepening of the electro- 
static potential profile occurs when the ion velocity Vi at the 
boundary exceeds the critical (Bohm) velocity vb ■ This leads 
to the acceleration of the negatively charged microparticles 
towards the center of the void. The electrostatic potentials 
on the void boundary before and after the critical transfor- 
mation are (l>vi and (I>v2 respectively. 

in an ideal case of a coUisionless sheath, corresponds to 
the Bohm velocity vb) in order to compensate it. Thus, 
our principal hypothesis is that the heartbeat instability 
occurs due to this critical phenomenon. 

As illustrated in Fig. |6j sudden steepening of the elec- 
trostatic potential profile occurring at the critical transi- 
tion has two important consequences. First, electric field 
at the void boundary increases dramatically, leading to 
the appearance of an uncompensated net force, which 
pushes the microparticles at the boundary towards the 
center. Second, in a transient situation of the sheath for- 
mation electron and ion fluxes through the void boundary 
are not necessarily equal. Electrons can therefore accu- 
mulate inside the void, causing the spatial redistribution 
of the glow intensity. Obviously, these consequences qual- 
itatively explain the observed dynamics of microparticles 
and intensity of the plasma glow inside the void (Fig.jsjj) 
during the heartbeat instability. 



B. Ionization balance inside the void 

Inside the void electrons and ions are produced by the 
electron impact ionization of neutral argon atoms, with 
the rate 

P 

Qi=n-kg,{Te), (1) 

where kgi is the rate constant of the electron impact ion- 
ization from ground state [33] and n is the number den- 
sity of bulk quasineutral plasma. The production rate of 
electrons and ions is compensated by their flux through 
the void boundary. Close to the critical condition of 
sheath formation the ionization rate should be equal to 



5 



the Bohm rate of ion loss, 



Qb = risVB 



(2) 




where = ne~^/'^ is the ion density at the sheath edge, 
Sy is the surface area of the void boundary, is the 
volume of the void and vb — \fTeJrni with mi, being 
the argon ion mass, is the Bohm velocity. This yields 



(3) 



where i?„ = V^jSy is the effective void radius. Equa- 
tion ([3]) determines the electron temperature, at which 
the sheath on the void boundary should be formed. Sub- 
stituting all values into Eq. ([s]) we obtain Tg = 3.1 eV 
for p = 20 Pa. 

The balance of ion and electron losses through the 
void boundary, UsVb = l/AnvT^e~'^f^^', where vt^ — 
yJWr\^Jwm^ with me, being the electron mass, is the ther- 
mal velocity of electrons, determines its floating potential 
(^/, which is equal to ~ 6.3Te for argon [21]. In a steady 
state, void boundary acquires the potential since the 
net dc current is zero. However, evolution of the glow in- 
tensity in Fig. [3] suggests that after the sheath is formed, 
transient currents should be induced in the system: void 
boundary should charge to a slightly larger negative po- 
tential (compared to ^y) in order to partly reject the 
electron flux from the void into the complex plasma. 

Let us assume, that the electron temperature remains 
constant upon the critical transformation at the void 
boundary. The observed increase of the emission inten- 
sity inside the void can then only be attributed to the 
variation of electron density, which in turn can be ex- 
pressed as follows: 



(Te) 



1 



exp 



(4) 



where Scjjf is the deviation of the potential of at the void 
boundary from (f>f and r is the rise time of the emission 
intensity in Fig. [3b. Also, 
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where 51 is the observed increase of the emission inten- 
sity I. After combining Eqs. Q and we notice that 
electron density cancels out and <5^/ can be expressed in 
the following way: 
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(6) 



For the value of = 3.1 eV, obtained above, Eq. ([6| 
yields 54)f/Tf. = 8 x 10~^. Therefore, even small devia- 
tions of the void boundary potential from the equilibrium 
(floating) value can cause the observed variations of the 
glow intensity. 



Number density of IS5 metastable state was found 
from the absorption measurements to be equal to nis^ = 
1.6 X 10^^ m^"^. We can define the ratio of the rates 
of electron impact ionization out of the metastable and 
ground states, 



s(T, 



p/T kgi{Te) 



(7) 



where kmi is the rate constant of electron impact ioniza- 
tion from the metastable state, expression for which can 
be taken from Ref. ^IT, For our case we obtain s ~ 10"''. 
Hence, ionization from metastable state is therefore con- 
sistently omitted in Eqs. ([T]), ([3| and Q. It plays, how- 
ever, the key role in the laser excitation of the heartbeat 
instability, making possible a small modulation of the 
ionization rate by the chopped laser beam (see Sec. IV E I . 



C. Dynamics of microparticles near the void 
boundary 

According to our model, it is the same critical trans- 
formation of the electrostatic potential profile at the void 
boundary, which is responsible for both the increase of 
the plasma emission inside the void and the acceleration 
of microparticles towards the center of the void. Then 
the floating potential </)/ can be related to the dynam- 
ics of microparticles. After the formation of the sheath 
microparticles will be accelerated by the potential differ- 
ence A^i, = \4>y2 — 4>vl\^ where 4)^2 and are the void 
potentials after and before the critical transformation, 
respectively (see Fig. [6]) , so that 



(8) 



Here m^;, v and are the mass, velocity and charge of 
a microparticle, respectively. The microparticle charge 
can be expressed as follows [?7] : 



(9) 



where z is a dimensionless factor depending on the ratio 
of electron and ion temperatures (and densities) as well 
as on the Mach number of streaming ions. We note that, 
according to the Bohm criterion [26j . the sheath forms 
when the drifting ions reach the energy of (I)b = 
Therefore, (pyi cannot significantly exceed the value of 
(pB- Then (p.^i (and, of course, 5(j)f) is substantially 
smaller than ipf ~ 6.3Te and, consequently, A^^, is some- 
what smaller than (j)f. Combining Eqs. ([8| and (joj) we 
obtain 
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Then, after the substitution of the corresponding param- 
eters in Eq. ( 10 ) and assuming z = \ deep inside the 
sheath we obtain v ~ 1.5 m/s. 



We measured the velocity of the accelerated micropar- 
ticles averaged over about 5 ms after they start their 
motion into the center of the void. Effect of the Epstein 
neutral drag [35], whose characteristic time for the con- 
ditions of our experiment is 13 ms, is then negligible. 
According to the measurements v ~ 0.5 m/s, which is 
not inconsistent with the value obtained from Eq. (10 1. 



D. Dependence on the discharge parameters 

We supposed that the heartbeat instability occurs due 
to the formation of the sheath at the void boundary. This 
happens when the diffusion rate of the plasma through 
the void boundary, Qdiff ~ nDa/Ry (where Da oc T^/p 
is the ambipolar diffusion coefficient), reaches the value 
of Qb determined by Eq. ^ . Due to the ionization bal- 
ance Qdi f f , should always be equal to Qi , determined by 
Eq. ([T]), which yields the following scaling: 



pRy cx jTe/kg,{Te). 



(11) 



We note that this should not be understood as a self- 
consistent scaling for the void size, since the balance of 
forces acting on the microparticles on the void boundary 
is not taken into account in Eq. (111. Next, we consider 
the ratio 



C 



Qdiff 

Qb 



CX cx 

Qb pRv 



kg.{Te), (12) 



which controls the formation of the sheath at the void 
boundary and, hence, the appearance of the heartbeat 

is a decreasing function 



instabiUty. The rhs. of Eq. (11 



of Te (at least while Tf, is less than the ionization poten- 
tial). This means, that the ratio C decreases with the 
product pRv 

Scaling (12) can be used to understand e.g. certain 



features of the stability diagram shown in Fig. [2j Ac- 
cording to Ref. [22] Tf. monotonically decreases with the 
discharge power in a microparticle-free rf discharge. This 
can explain the existence of the upper power threshold 
for heartbeat instability. On the other hand, monotonic 
decrease of Te with the discharge power contradicts the 
existence of the lower power threshold. However, in a dis- 
charge with microparticles this dependence could easily 
become non-monotonic due to e.g. variable microparticle 
charge or size of the void. Furthermore, since the param- 
eter C decreases with pRy , presence of a pressure thresh- 
old in Fig. [2] becomes also clear. Increase of p suppresses 
the ambipolar diffusion, so that above some critical value 
of pressure the sheath on the void boundary can be no 
longer sustained irrespective of the rf power. 

In the framework of this approach, the onset of the 
instability upon the increasing of the number of levitat- 
ing particles can also be explained: Injecting additional 
microparticles into the discharge leads to the decrease of 
the void size Ry and, therefore, increase of the ratio C. 
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FIG. 7: Fourier spectra of the light intensity inside the void, 
for the experimental points from Fig. |4^. At fias = 46 Hz 
(when the intensity of the heartbeat is close to maximum) 
the spectrum is peaked at about the half of the laser fre- 
quency, suggesting the parametric instability as underlying 
mechanism of the laser-excited heartbeat. 



E. Effect of the tunable laser beam 

The laser-excited instability is not the result of the 
radiation pressure exerted on the microparticles. First 
of all, the instability was observed at the laser intensity 
of J = 600 W/m^, which is hardly enough to produce 
significant dynamical effects [SO]- Also, the instability 
disappears once the laser wavelength is detuned from the 
center of the absorption line. The effect produced by 
the laser is, therefore, due to its interaction with the 
background plasma. 

The fact that the instability can only be triggered 
when the laser beam hits the void indicates the minor 
importance of any kind of direct interaction between 
metastable atoms and microparticles (e.g., secondary 
electron emission [3T]). Therefore, in case of laser ex- 
citation it is only the perturbation of the ionization that 
can lead to the instability. Ionization from metastable 
states provides only a small fraction s ^ 10^^ to the 
net ionization rate. This suggests that the observed res- 
onance is the consequence of the parametric instability. 
To address this issue we analyzed the Fourier spectra 
of signals from the photomultiplyer, pointed to the void 
(Fig. [7]). Close to the resonance at fias = 46 Hz, when 
intensity of the heartbeat is high, we observed excitation 
of a broad frequency band with the peak in the vicin- 
ity of 23 Hz, i.e., halved driving frequency. This is one 
of the characteristic features of a parametric instability. 
At fias = 49 Hz, when no heartbeat is observed, this 
band is very much suppressed, as compared to the case 
of fias — 46 Hz. Due to the parametric instability the 
plasma parameters inside the void could be changed in 
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such a way, that the sheath formation at the void bound- 
ary could be triggered. 

Next, we show that the laser beam with the intensity J 
tuned to the wavelength of 2p^ — >■ Iss transition of argon 
can significantly disturb the balance of metastable atoms. 
The laser light is absorbed by the Iss argon metastable 
atoms is absorbed with the rate 

Qabs — r^abs^ilss- (13) 

Here F = XJ/hc is the flux of laser photons and 



8^ 



92p- 



2T gu 



A. 



(14) 



is the absorption cross-section, where, g2p^ = 3 and 
gis^ = 5 are the statistical weights of the respective 
atomic states and A = 5.18 x 10^ s~^ is the Einstein 
coefficient of the transition 132 



Equation 14 holds only 
for the center of a spectral line, with a Doppler profile 
characterized by the temperature T. Absorption mea- 
surements, performed with the same laser and using the 
same transition, showed that the spectral line indeed 
has a Doppler shape with T « 300 K. By substituting 



these numbers into Eqs. (13) and (14), we get Qabs = 
5.6 X 10^"'^ m~^s~^. Branching fractions for the atomic 
level 2p7 of argon are such that only 23% of atoms, ex- 
cited to this level return into metastable states [35] . The 
rest will decay to ls2 and IS4 radiative states, whose life- 
times are orders of magnitude shorter than those of the 
metastable states. Therefore, our laser beam produces 
an effective sink of metastable atoms. 

For comparison we estimate the rate of electron impact 
excitation from ground state into IS5 metastable state for 
the typical value of electron density Tig — 10^'* m~'^ and 
temperature — 3.1 eV (obtained from the ionization 
balance). The rate of this process, Qexc = 3.4 x 10^^ 
m~''s~^ [23] is much smaller than the rate of destruction 
of metastable states by the laser beam. Number density 



of metastable states inside the laser beam is significantly 
reduced compared to the unperturbed plasma. 



V. CONCLUSION 

The heartbeat instability is identified as a heteroge- 
neous phenomenon, occurring when a clear interface ex- 
ists between the microparticle-free region (void) and com- 
plex plasma. In the proposed simple model, the instabil- 
ity occurs due to the formation of the sheath on the void 
boundary. Abrupt change of the electrostatic potential 
profile leads to the appearance of uncompensated electro- 
static force pushing the microparticles towards the center 
of the void. Transient variations of the potential of the 
void boundary can be responsible for the observed varia- 
tions of the plasma glow intensity in the entire discharge. 
We have also shown that our hypothesis can be useful 
in understanding the effect of different parameters {i.e. 
discharge power, gas pressure and number of levitating 
microparticles) on the heartbeat instability. Excitation 
of the heartbeat by a modulated tunable laser beam was 
found to be the result of a parametric instability. 

To conclude, we would like to stress that our hypothe- 
sis on the mechanism of the heartbeat instability should 
not be considered as a self-consistent approach, aiming to 
completely solve the problem. Clearly, careful dedicated 
experiments are needed to further confirm or deny this 
hypothesis. 
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